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• We present the diagrammatic technique for calculating the free energy of the matrix eigenvalue 
I model (the model with arbitrary power /3 by the Vandermonde) to all orders of 1 /N expansion in 

the case where the limiting eigenvalue distribution spans arbitrary (but fixed) number of disjoint 
intervals (curves) and when logarithmic terms are present. This diagrammatic technique is corrected 
and refined as compared to our first paper with B.Eynard of year 2006. 

> 

o , 

1 Introduction 

ON 

The standard Hermitian one-matrix model is determined, after integrating out angular degrees of 

• freedom, by the A^-fold integral over eigenvalues of the form 

o ■ 

N 



J fl dxi A(x)2e-^ h = to/N 

1=1 



where A(x) is the Vandermonde determinant of eigenvalues Xi, h is the formal expansion parameter, 
and to is the normalized number of eigenvalues. The /3-ensembles, or /3-eigenvalue models, are gener- 
alizations of this integral obtained by replacing the exponent of the Vandermonde determinant by an 
arbitrary positive number 2/3. In the proper normalization, we then evaluate the (formal) integral 

N 

ncix,|A(x)|2/3e-^^£i^(-') = e-^. (1) 

i=l 

The three values /? = 1/2, 1, 2 correspond in the Wigner classification by the respective ensembles of 
orthogonal, Hermitian, and self-dual quaternionic matrices, but we can put forward the problem of 
calculating the perturbation expansion of the /3-model free energy for arbitrary value of (3. 

In our joint paper with B. Eynard [5], we developed the perturbation expansion solution to the 
integral ([T]). Our procedure enables constructing all the correlation function and the free energy 
expansion terms order by order of the double asymptotic expansion in the two parameters: h"^ and 
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{VP - The terms of these expansions (for all the terms except few leading terms) were 

expressed using the special Feynman-like diagrammatic technique, which is a generalization of the 
diagrammatic technique originally constructed in our paper [4J for the Hermitian one-matrix model. 

The recent interest to the /3-eigenvalue models is due to the Alday, Gaiotto and Tachikawa conjec- 
ture [1] relating instantonic generating functions to the conformal blocks of the Liouville theory (see 
also [9] ) . These conformal blocks were in turn associated with the asymptotic limits of the /S-ensembles 
containing logarithmic terms in the potential. More precisely, it was shown (see [7] and references 
therein) that the Penner-like model with three logarithmic terms proposed by Dijkgraaf and Vafa [6] 
correctly reproduces the standard results of = 2 supersymmetric gauge theories. 

On the other hand, explicit calculations performed by Sara Pasquetti and Marcos Marifio had 
revealed that applying the original diagrammatic technique of [5] results in a nonsymmetric result for 
the two-point correlation function already in the subleading expansion order. This clearly indicated 
the incompleteness of the diagrammatic technique of our first paper. 

Therefore, the aim of this paper is twofold: on the one hand, we extend the pattern of the paper [5] 
including the logarithmic potentials into consideration and, on the other hand, we formulate a new 
improved diagrammatic technique endowed with the new selection rules that excludes part of diagrams 
that caused a mismatch when calculating correlation functions. 

As a nearest perspective, it is interesting to include the hard edges into consideration, i.e., to 
construct the consistent diagrammatic technique for the /3-ensembles in the situation where we put 
explicit restrictions on the intervals of the eigenvalue distribution; the first step in this direction was 
done in the recent paper [2] in which a one-cut solution was investigated and results for few lower 
order terms of the free-energy expansion were derived. 



2 Eigenvalue models in expansion 

In this section we show how the technique of Feynman graph expansion elaborated in the case of 
Hermitian one-matrix model [8] , [1] can be adapted to solving the (formal) eigenvalue model with the 
action 

„ TV t / N \ 

i=l k=l \i=l J 

where V{x) = Y^s>i^sX'^ , ^ = to/^ is a formal expansion parameter, and tQ is the normalized 
eigenvalue number. The integration in ([2]) goes over N variables Xi having the sense of eigenvalues 
of the Hermitian matrices for /3 = 1, orthogonal matrices for /3 = 1/2, and symplectic matrices for 
f3 = 2. In what follows, we set /3 to be arbitrary positive number. The integration may go over curves 
in the complex plane of each of A^ variables Xj. For /3 7^ 1, no topological expansion in even powers 
of h exists and we rather have the expansion in all integer powers of h. We also assume the potential 
V{p) to be a polynomial of the fixed degree m, + l, or, if we include logarithmic terms in the potential 
introducing 

t 

y(xi) = y(xi) ^7felog|xi - Ofcl, (3) 
fe=i 

then it suffices to demand the derivative V'{p) to be a rational function [8]. 

For brevity, we let the symbol tr/(A) to denote X^i^i fi^i) ^^i any function f{x). 
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2.1 The loop equation and resolvents 

The averages corresponding to partition function ([2|) are defined in a standard way: 



(HX)) = \ I^DXFiX) exp (^-^trF(X)) 



(4) 



with the total potential V{x), and we introduce their formal generating functionals: the one-point 
resolvent 

A;=0 ^ 

as well as the s-point resolvents {s > 2) 



k„tks = l Pi ■■■Ps \ Pl-X PS-X/ 

where the subscript "conn" pertains to the connected part. 

These resolvents are obtained from the free energy through the action 



W{pi,...,ps) = -h 



2 d d dF 



dV{ps) dV{ps-i) dV{pi) 
' ' ' Wipi), (7) 



dvips) dv{ps^i) dv{p2) 

of the loop insertion operator 

dV{p) " ^ pJ+i dtj • ^ ' 

Therefore, if one knows exactly the one-point resolvent for arbitrary potential, all multi-point resolvents 
can be calculated by induction. 

The loop insertion operator has an interpretation in the theory of symmetric forms. Let us set into 
the correspondence to an s-point correlation function Wp^^,,,^p^ the symmetric form • • • dps 

that is meromorphic on our Riemann surface. We let Of q denote the (infinite-dimensional) space of 
such forms. Then, the action of the loop insertion operator increases the order of this form by one, 

A : ^ (9) 

We now consider the perturbation expansion of quantities in ([7]). In the above normalization, the 
/i-expansion has the form 

oo 

W{pi,...,ps) = Y.K'Wri2{pi,...,Ps), s>l, (10) 

r=0 

and for /3 7^ 1 it comprises both even and odd powers of ^, whereas in the Hermitian matrix model 
case (/3 = 1) the corresponding expansion (the so-called genus expansion) goes over only even powers 
of h. 

We obtain the master loop equation of the eigenvalue model ([2]) if we perform the changing Xi — >■ 
Xi + ^ of the integration variable. The resulting (exact) equation first takes the following form in 
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terms of the means 

N 



i'j^^\(.Xi-p){xj -p) / f^^\{xi-p) 

N /m+1 



i=l \ s=l 



Xi-p 



^^^/^ 1 

' n mi. — n 



k=l 1=1 



i=l 
1 



Xi-p 



Xi Oij; 



Xi-p 



0. (11) 



We let Pm-iip) denote the degree-(m — 1) polynomial X^i^i (^127=1 ■s^s ^' 3.-„p — ^ Note also that the 
resolvent W{p) must be nonsingular at p = for any k. Then Eq. pip becomes in normalization 
(©: 



(12) 



w\p) - h'w{p,p) + ^ V/3 - V/3 w'{p) + y'(p)w(p) + Pm-iip) 



* 1 

V7fc [W{ak) - W{p)] = 0. 

nji. — r) 



fc=l 



OLk-P 



2.2 Solution at large N 

Here we specially consider the limit as — in (fT2]l . Note that VF^°-'(p)|p^oo = ta/p-'r 0{p~^). 

For VF('^^(p) we have the algebraic equation to resolve which it is useful to introduce the function 



y{p) = W{p) 



v'{p) - 



ik- 



k=l 



P-Oik 



W{p) - -V'{p) 



(13) 



for which Eq. (jl2p gives the full square: 



y{p) 



^/P2m+2t{p) 

nLib-«fe) 



(14) 



The function y{p) may have up to 2m + 2t branching points thus defining a hyperelliptic Riemann 
surface of maximum genus m + t — 1. We let 2n be the actual number of the distinct branching 
points fir (some branching points may merge producing double points), and the actual genus is then 
n — 1. We assume that we have therefore n disjoint intervals Ai of eigenvalue distribution and n — 1 
"forbidden zones" Bi . We treat n — 1 among n intervals Ai as a-cycles and we then let X] j=/ 
be the corresponding 6j-cycle (one half of which goes along the physical sheet and the other half goes 
along the second, unphysical sheet). So, in general, we assume that y(p) has the form 



y{p) = Mm+t^ 



-n{p) 



Mr 



nLib-afc) 



where the polynomial Mm+t-n{p) arises from merging branching points of P2m+2t{p)- 

The ra + t — n + 1 coefficients of the polynomial Mm+t~n{p) and the 2n branching points are 
to be determined from the following conditions: 



(asymptotic condition at infinity) 



i 



^lim y(p) = -iy'(p)+ |^^2^7fc + io 



P 



+ 0{p' 



m + 2 conditions; 
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(regularity at p = of W{p)): 

■\/ -f2m+2t(aA:) _ 1 



7fe t conditions; 



(equality of chemical potentials on all cuts) 

ydj) = n — 1 conditions. 



/ 

JB, 



Note that, in the context of the Seiberg-Witten theory, the last n — 1 conditions can acquire other 
forms, for instance, 

• (fixing filling fractions on the intervals of eigenvalue distribution) 

ydp = €i n — 1 conditions. 



/ 



So, in total wc have exactly the desired number ofm + t + n + 1 conditions. 

The recurrent procedure in the subsequent sections is such that all the geometry properties of 
the theory (Bergmann kernels, etc.) is completely determined by the reduced hyperelliptic Riemann 
surface 



2n 



r=l 

where we lose all the information on the rational function Q{p) := Mm+t-n{p) / iYl\=i{p — ctfc))- 

All the rest of the text is totally insensitive to the nature of the rational function Q{p). 

Wc need to introduce the Bergman kernel: the symmetric {B{p,q) = B{q,p)) bi-differential that 
has the form 

dxi'Di dx { Oi f 

B{p,q) := B{x{p),x{q))dx{p)dx{q) —— ---^ as x{p) ^ x{q), f B{-,p)=0 

{x{p) - x{q)) 

in the local coordinate x(j)). 

We also need its primitive (a one-differential in r): 



dEp,q{r):= / i?(-,r) 

Jx(n) 



x(p) 

Note that both B(p,q) and dEp_q{r) depend only on "reduced" Riemann surface (are determined 
completely by the branching points fir)- Definitely, they take the simplest form for one-cut solutions: 
if we have just two branching points /xi and iJ,2, then, for q being the image of q on the other sheet. 



dE -(r) - 

for higher-genus spectral curves, this expression contains also projections terms proportional to Abelian 
differentials normalized to i?-cycles. 
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2.3 Higher-order corrections 

The /3-dependence enters (jl2p only through the combination 



C=7^-V^, (15) 
and, assuming /3 ~ 0{1), we have the free energy expansion of the form 

oo oo 

J- ^ T{h, C, to, ti, i2, . . .) = E E ^''^'''C'^k,i. (16) 

k=0 1=0 

Substituting expansion pO]) in Eq. p^ . we find that Wg(p) for 5 > 1/2 satisfy the equation 

dV{p) dp 



g'=l/2 



+P^Up) + i2^k^—W,ia,) (17) 



In Eq. (flTl) . Wg(p) is expressed through only the VF^. (p) for which gi < g. This fact permits developing 
the iterative procedure. 

In analogy with ([16]), it is convenient to expand multiresolvents Wg{-) in Q: 



Then, obviously, (fTTjl becomes 



-2y{p)Wk,i{p) = E ^fci,;i(^^)^^-fci,/-h(p) + ^^^fe-i,Kp) + C^M^M-i(^') 



*;i + !l>0 



Pin%) + il^s-^WuA0Cs). (19) 



Recall that because 



.=1 ^'-"^ 



VFfc,/(p)|p^oo = -5fc,o<5/,o + 0(l/p2), (20) 



all Wk,i{p) are total derivatives. 



^^M=g^fk,U k,l>0. (21) 



The planar limit solution for VFo,o(p) exactly coincides with the one in the one- matrix model. The 
normalizing conditions seem to be those ensuring the coincidence of chemical potentials in all orders 
of the perturbation theory 



f Wk,i{C,P2,---,Ps) = 0, for k + l>0. 



(22) 



Note that in the Seiberg-Witten theory, these conditions must be merely replaced by the vanishing 
conditions for the j4-cycle integrals. 
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The formal solution to ()19p is given by the integral [5]: 

''^D yy^J 0<ki<k,0<li<l 

0<ki+li<k+l 

Here the contour lies on the physical leaf, encircles all cuts of the function y{q), and leaves outside 
all the other possible singularities (poles and zeros). In the next section we develop the diagrammatic 
representation for terms of these recurrent relations. 



3 Graphical representation. Spatial derivative. 



The diagrammatic technique of [5] was based on the diagrammatic technique developed in [1] for the 
Hermitian one-matrix model, but contained new elements. The original technique of [3] contained two 
types of the propagators: arrowed lines corresponding to dEq^q{p) and nonarrowed lines corresponding 
to B{p,q). Simultaneously, the two-point correlation function 



d 



dV{p) 



Wo,o{Q) = Bip,q) 



1 



(p - 



and the action of the loop insertion operator on it produces the three-point correlation function 
described by the Ranch variation formulas that result in the appearance of the three-point vertex 
denoted by light circle and containing the integration over the contour as in ()23p with the weight 
l/(2y(g)) First of all, the recursion kernel in ()23p coincides with the recursion kernel in but to 
comply with the new term, C^^fe,«-i(p)j appearing in the right-hand side, we must introduce a new 
propagator dpdy{q) denoted by the dashed line. 

The graphical representation for a solution of loop equation (jl7p then looks as follows. The multi- 
point resolvent Wk^i{pi-, . . . ,ps) is represented by the block with s external legs and with the double 
index k, I. Among these legs one, say pi, is selected to be the root of maximal tree subgraph composed 
from arrowed propagators and establishing a partial ordering on the set of vertices of a graph. 

We present the derivative -^Wk,i-i{pi, ■ ■ ■ ,Ps) as the block with s + 1 external legs, one of which 
is the dashed leg that starts also at the vertex pi. For instance, for the one-point resolvent expansion 
term Wk^i{p), we obtain [5]: 



p Q 



E 

^1 




p Q 




(24) 

In the case of multipoint resolvent Wk^i{pi, ■ ■ ■ ,Ps) the picture is similar, we replace p by pi and add 
s — 1 remaining external legs to the diagrams in the right-hand side; then, in the first diagram we 
must distribute the added s — 1 external legs over two subdiagrams in an arbitrary way and take a 
sum over all possible insertions that leave the subdiagrams Wk^^i^iq^Pj^, ■ ■ ■ ^Pjs^) "stable." Recall that 
Wkiipii ■ ■ ■ ,Ps) is called stable (for a nonempty set of pi, . . . ,ps) if either /c > or / > or s > 2. 

We can then in turn present the term Wk^i~i{q,P2, ■ ■ ■ ,Ps) in the form ()23p with the recursion 
kernel dErj^fjiq) / C^uiv)) integrated with some function F{rj), the rest of the diagram, which we do not 
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specify here. Let us consider the action of the spatial derivative d/dq on Wk,i-i{q-,P2-, ■ ■ ■ ,Ps)- We 
place the starting point of the dashed directed edge at the point q and associate just dx with this 
starting point. The first object (on which the derivative actually acts) is dErj^ftil), then comes the 
vertex with the integration ^ §Av) "t^; then the function F{ri). We can present the action of the 
derivative via the contour integral around q with the kernel B{q,(^): 

B_(t dE.,„-,h)d, \ / Bi,.i),E,.,ii) 

dq\Jc'n) 2^iy{,,) ^"J (^,JcM d( 2iri 9(7,) ^ ^ ' 

where q lies outside the integration contour for r]. The integral over ^ is nonsingular at infinity, so we 
can deform the integration contour from Cq to > 

B{q,C)dE,,^{0 j; d^ -F(r?) = -/ B{p,,OdE^^^[0 dy ^^^^ ^^^^ 



2Tri d^ JcM2my{ri) Jc^^^>d^^ 2m d^ 27ri y{r])' 

We now push the contour for ^ through the contour for ij, picking residues at the poles in ^, at ^ = rj, 
i = f] and at the branch points. We then obtain 

^^^B(,.()iE,,(() .(^) J 3S4-Fi,,)+l ^?M>F(,), (27) 

M^£^|.. <i{ 2Ti >,(,,) ^ " /(.M 2»i !/(,,) " JcM> 2m y(r,) ^ ^ ' 

where the first integral has only simple poles at the branch points. The main trick stems from that the 
residue remains unchanged if we apply the I'Hopital rule and replace B{q,Q / dy{S) by dE^^^{q)/2y{^), 
thus replacing B{q,^) by dE^^^{q)dy{^)/2y{^^ we then obtain in the right-hand side of 

^p^,, dE^,^^P^)dy{OdE,,,{0 d^ r B{q,^)-B{q,f^) 

W^C^M. 2y{0dC 27riyirj) ^'^ fcM 2m y{i^) ^'^^ 



Pushing contour of integration for around the branch points back through the contour for rj, we pick 
residues at ^ = and £, = fj, which both give the same contribution, and obtain 

dE^^^{q)dy{OdE^rj,f,{0 drj ^ 
c(^)^C^v) 2TTi2y{Od^ 27riy(r?) ^"^^ 



W y{rf)d'q 2Tiiy{ri) Jc^n) 2iTi y{rj) 

We evaluate the last two terms by parts therefore obtaining 

dE^^^{pi)dyiOdE^,f^iO dr^ ^ r dE,,fiiq)drj dFjy) 

ci^Hc^^^ 2m2y{i)di 2Tri y{ri) ^'^^ JcM 2m y{ri) Or] ' 



We have thus found that 

d I f dEr^^n{q)dri 



9q \ Jc^i!]^ 2m y{7]) 



dE^-^{q)y'{i)dE^^^{i) dy . , / dE,,f,{p,) dr, dFjy) 

c(«)>c(,'') 2m2y{0 2m y{ri) ^'^^ JcM 2m y{ri) Ot] ' 



'^Here and hereafter, the comparison of contours indicates their ordering. 
^We therefore "imitate" pole terms by inserting dy/y. 
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and we can graphically present the action of the derivative as follows: 

d.^olFiv)} = ^Wo{i^(r?)| + ^0{|F(r?)| 

q rj^ J q t T]^ > Q rj ^ ^ > 

(31) 

In these relations and in what follows, we use light circles ("white" vertices) to denote the integrations 
with the weight l/(2y(r/)) and "black" vertices to indicate integrations with weights containing the 
first or higher derivatives y^^\ri) in the numerator. 

We also assume the ordering from left to right in all the diagrams below. 

We therefore see that using relation ()3ip we can push the differentiation along the arrowed prop- 
agators of a graph. It remains to determine the action of the derivative on internal nonarrowed 
propagators. But for these propagators (since it has two ends), having the term with derivative from 
the one side, we necessarily come also to the term with the derivative from the other side if we act by 
the dashed line that was commenced before the outer end of this propagator. Combining these terms, 
we obtain 

and we can deform this contour to the sum of contours only around the branching points (sum of 
residues). Then we can again introduce y' {S,)dS,/y{C) and integrate out one of the Bergman kernels 
(the one that is adjacent to the point q ii p > q or p if q > p; recall that, by condition, the points p 
and q must be comparable). 

That is, we have 

r p q 




We must therefore improve the diagrammatic technique of /3-model in comparison with the Her- 
mitian one-matrix model by including the dashed lines. They are not actual propagators, but these 
lines ensure the proper combinatorics of diagrams. Indeed, from (jSip it follows that the derivative 
action on the "beginning" of the dashed line is null, dpdp = 0, and when this derivative acts on the 
"end" of this line with the black vertex, we merely have 

d,y'{q)=y"{q), (33) 

which we denote symbolically as two dashed propagators ending at the same black vertex. If we 
continue to act by derivatives d/dq, then, obviously, when k dashed propagators are terminated at 
the same vertex, we have the kth order derivative y^^\q) corresponding to them. 

Another lesson from studying the action of the spatial derivative is that not all diagrams are 
possible: for example, the diagram containing the vertex 




r p q r) 
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never appear in our approach. Later on we formulate all the selection rules for the diagrams as well 
as present all types of possible vertices. 

We also introduce the black-and-white coloring of vertices in accordance with the rule: if there is 
no factors y^^\q) standing by the vertex, it is white; if there are such factors, the vertex is painted 
black. 

We need now to produce the action of the loop insertion operator in the new setting. We calculate 
the action of d/dV{r) on the dashed propagator. Obviously, 

r) r)^ 

y^''Hq) = ^B{r,q), (34) 



and any attempt to simplify this expression or to reduce it to a combination of previously introduced 
diagrammatic elements fails. This means that we must consider it a new element of the diagrammatic 
technique. We indicate it by preserving k dashed arrows still terminating at the vertex q with the 
added nonarrowed solid line corresponding to the new propagator B{r,q). The vertex then change 
the coloring from black to white because it contains y^^\q) factors no more and we therefore assume 
that these derivatives act on B{r, q). 

An accurate analysis demonstrates nevertheless that we must include also vertices with the mixed 
action of spatial derivatives: on B{r,q) and on y^^\q) standing at this vertex (which is then also 
painted black). These vertices were also missed in the original paper |5j. 

We leave the detailed proof for subsequent publications and present in the next section the complete 
set of rules for constructing diagrammatic representation for Wk,i{pi, ■ ■ ■ ,Ps for all k and I. 



4 Feynman diagrammatic rules 



We now describe the diagrammatic technique that results from the action of dq — the spatial derivative 
and the action of the loop insertion operator d/dV{q) with describing new vertices and selection rules. 

In this diagrammatic technique, all the solid arrowed propagators dEq^q{p) are free of spatial deriva- 
tives, vertices always contain the factor l/(2y(g)) and may or may not contain additional derivatives 
y^^\q) in the numerator and the solid nonarrowed propagators corresponding to S-kernels may con- 
tain additional spatial derivatives acting only on the innermost end of the propagator (recall that 
-B-propagators can join only vertices comparable in the sense of partial ordering established by the 
tree of arrowed propagators. 

In the graphical representation below the ordering of vertices is implied from left to right. 

The diagrammatic technique for constructing resolvents Wk^i{pi, ■ ■ ■ ,Ps) is as follows. 

We take the sum over connected graphs containing k internal nonarrowed propagators (Bergmann 
kernels), / dashed arrowed propagators (those are always internal "derivatives") and s external legs 
such that 

• we begin constructing each graph by constructing the maximum rooted directed tree subgraph 
composed from the propagators dEq^q[p) [q here is the inner vertex with respect to p, the arrow 
is pointed toward q) and take the sum over all choices of such rooted tree subgraphs; 

• we choose the external vertex pi to be the root of the tree, one and the same pi for all the 
graphs; 

• all other external legs are nonarrowed propagators B] we call propagators dE and B solid as 
we set solid lines corresponding to them; all the external points of the graph are considered 
outermost points; 
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• other internal edges (i.e., edges that are incident only to the graph vertices and not to pi, . . . ,ps) 
are either B (there are exactly k of them) or dashed arrowed edges (there are exactly I of them) ; 

• rooted tree subgraph establishes a partial ordering on the set of internal vertices; we allow 
internal B edges and dashed edges to connect only comparable vertices. Dashed arrows are 
directed always inward, but can begin and terminate at the same vertex; if such a vertex exist, 
it must be an innermost vertex of the graph. The i?-lines can also begin and terminate at the 
same vertex; in this case, this vertex is also necessarily an innermost vertex of the graph. 

• A partial ordering establishes the order of taking residues: we start from innermost vertices and 
continues toward the root of the tree. 

• All s external legs pj are outside all the integrations. 

• All the internal vertices are incident to exactly one incoming arrowed line and can be incident 
up to two other ends of solid lines. 

• The vertices of the graph are painted either black or white (depending on whether they re- 
spectively contain or do not contain the factors y^''\q)) and must satisfy the selection rules 
below. 

• All the graphs satisfying all the above properties enter the final answer with the standard sym- 
metry coefficients that are the reciprocal volumes of the discrete automorphism groups of the 
graphs. 



The selection rules follow from the complete list of ten types of vertices below. Their origin will 
be partly clarified in the next section when considering the action of the loop insertion operator on 
the elements of the diagrammatic technique. The general proof will be published elsewhere. 

The vertices with three adjacent solid lines 




Here the vertex p must be outside all the starting points of dashed lines acting on the vertex q and 
can be external. 

(37) 
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S dE„ A* 



dq 



2m y{q) 



B{q,*)dE,4q) 



(38) 



(39) 



The vertices with two adjacent solid lines 



dq dP 
2'Ki y{q) dqP 



{B{r,q)y^'\q)), 



k>0, p>0. 

Here, note that k must be greater than zero, so the point r cannot be external, 



(40) 



dq 



2m y{q) 



y^''\q)dE^4q), k > 0. 



(41) 




(42) 

Here the vertex p must be outside all the starting points of dashed lines acting on the vertex q and 
can be external. 



/c^?' 27rz y{q) 



(43) 



The vertex with one adjacent solid line 
,k 



dq 



2m y{q) 



y('=+i)(g). A; > 0. 



(44) 



When calculating Wk^i{pi-, . . . ,Ps) the order of integration contours (previously the order of taking 
residues at the branching points in the case of the Hermitian one- and two-matrix models) is prescribed 
by the ordering of vertices in the tree subgraph: the closer is the vertex to the root, the more outer 
is the integration contour. In contrast to the Hermitian matrix model case, the integration cannot 
be reduced to taking residues at the branch points only; all internal integrations can be nevertheless 
reduced to sums of residues, but these sums may now include residues at zeros of the additional 
polynomial M{p) on the nonphysical sheet and, possibly, at the point oo_. 

As in the Hermitian matrix-model case, we use the -ff-operator constructed below to invert the 
action of the loop insertion operator and obtain the expression for the free energy itself. 
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Figure 1: Cuts in the A-, or "eigenvalue," plane for the planar limit of the /3-model (coinciding with 
the planar limit of the Hermitian one-matrix model). The eigenvalues are supposed to be located "on" 
the cuts (curly lines). We add the logarithmic cut between two copies of the infinity on two sheets of 
the hyperelliptic Riemann surface in order to calculate the derivative w.r.t. the variable to; Cl is the 
corresponding regularized integration contour. 



5 Inverting the loop insertion operator. Free energy 

The ff-operator. We now introduce the operator (the functional on 1-forms) that is inverse to 
loop insertion operator ([8]). For a form df{x), leiH 

H^df{x) = -veSoo+V {x)df{x) - -veSoo^V {x)df{x) - to / df{x) -Y^Sii df{x). (45) 

The arrangement of the integration contours is as in Fig. [TJ 

The i/-operator, in contrast to the loop insertion operator, therefore reduce by one the degree of 
the (symmetrized) form associated with an s-point correlation function 

H : QI q ^ (46) 

We calculate the action of H on the Bergmann bidifferential B{x, q) using again the Riemann bilinear 
identities. We first note that as B{x,q) = dxdEx^qf^{q), we can evaluate residues at infinities by parts. 
Then, since dEx^qg^q) is regular at infinities, we substitute 2y{x) + 2to/x for V'{x) as x ^ oo+ and 
—2y{x) + 21^1 X for V'{x) as x — )• oo_ thus obtaining 

-reSoo+ f y(x) + ^ j d^^,g(,(g)da; + reSoo_ \ -y{x) ^^-^ dEx^q^{(i)dx 



_ n— 1 „ 

-todE,,q^{q)''~"^ -y2Si(f B{q,x). (47) 

i = l ''^^ 

Whereas the cancelation of terms containing is obvious, it remains only to take the combination 
of residues at infinities involving y{x). For this, we cut the surface along A- and -B-cycles taking 
into account the residue sX x = q. The boundary integrals on two sides of the cut at Bi then differ 
by dEx^qoiq) — dE^^^^ g^^{q) = 0, while the integrals on the two sides of the cut at Ai differ by 

dEx^qo{q) — dE^^^^ q^{q) = <f^ B{q,x), and the boundary term therefore becomes 

n-1 



i y{x)dx (b B{q,£,), 



which exactly cancels the last term in (j47p . Only the contribution from the pole at x = g then survives, 
and this contribution is just —y{q). We have therefore proved that 

H.,-Bix,q) = -y{q)dq. (48) 



^This definition works well when acting on 1-forms regular at infinities. Otherwise (say, in the case of VFo(p)), the 
integral in the third term must be regularized by replacing it by the integral along the contour Cl depicted in Fig. [T] 
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Let us now consider the action of the operator given by (j45p on Wk,i{x) subsequently evaluating 
the action of loop insertion operator on the result. Due to the commutation relations between and 
d/dV{x), 

H o d/dV — d/dV o H = 1A modulo the kernels, 

we have, 

{H^ ■ Wk,i{x)) = Wk,i{p) + ■ Wk,i{x,p)- (49) 

For the second term, due to the symmetry of Wk^i{p,q)-, we may choose the point p to be the root 
of the tree subgraphs. Then, the operator Hx always acts on (or, possibly, on its derivatives 

w.r.t. ^) where ^ are some integration variables of internal vertices. 

Let us consider the action of d/dV{p) on the elements of the Feynman diagram technique in SeclH 
Here we have three different cases; in all of them the variable p is external. 

• When acting on the arrowed propagator followed by a (white or black) vertex, we use the 
standard variational relation in ()50p below, where, again, -F(ry) is any combination of elements 
of the diagrammatic technique, 



d 



0{F(r?)} = ^^o{i^(r/)} 



dV{p) q r]^ J q ^ r] 



When acting on nonarrowed internal propagator ■^B{r, q), r > q, k > 0, we apply the relation 



(50) 




d 

dV{p) dq^ T^^^^q 

- - (51) 
without subsequent expanding the action of the derivative ^ into a sum of diagrams. 

Eventually, when acting on the factors y^^\q) standing at the black vertices, using relation (|34p 
we obtain from the diagrams (l40l) . (j4T]l . and (f44]) the respective diagrams (135]) . (f36]l . and (f42l) . 



We now consider the inverse action of the ff-operator in all three cases. 

In the first case where it exists an outgoing arrowed propagator dEfj^fj{(^) we can push the integration 
contour for ^ through the one for p; the only contribution comes from the pole at ^ = p (with the 
opposite sign due to the choice of contour directions in Fig. [1]). Graphically, we have 



^^^^^o{i^(??)| = ^o{f{v)} 

7 £ r/ J q 1]^ J 

(52) 



In the second case, the vertex ^ in ()5ip is an innermost vertex (i.e., there is no arrowed edges 
coming out of it). The 1-form y{£,)d(, arising under the action of H (I48p cancels the corresponding 
form in the integration expression, the expression becomes regular at the branching point and the 
residue vanishes. Graphically, we have for k >0, 



0. 

(53) 
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Eventually, in the third case, the action of the ff-operator just erases the new i?-propagator with 
the external variable p and restore the term y^''\q), and from the diagrams (j35p . (j36p . and (j42p we 
return to the respective initial diagrams (j40p . (j4ip . and (j44p changing back the color of vertices from 
white to black. Note that in this case we have the plus sign upon the action of the operator H. 

For Hx ■ Wk,i{x,p) = Hx ■ av^(^) Wk,iip), we obtain that for each solid arrowed edge, on which the 
action of d/dV{x) produces the new (white) vertex, the inverse action of Hx gives the factor —1, on 
each nonarrowed edge, on which the action of d/dV{x) produces the new vertex accordingly to (|5ip . 
the inverse action of Hx just gives zero, and at each black vertex, at which the action of the loop 
insertion operator changes the color to white and adds a new i?-propagator, the inverse action of Hx 
gives the factor +1. 

As the total number of arrowed edges coincides with the total number of vertices and the contri- 
butions of black vertices are opposite to the contributions of arrowed edges, the total factor on which 
the diagram is multiplied is minus the number of white vertices, which is exactly 2k + I — 1 for any 
graph contributing to Wk^i{p)- (For the s-point correlation function Wk^i{pi, ■ ■ ■ ,Ps) this number is 
2k + l + s - 2.) We then have 

Hx ■ Wk,i{x,p) = -{2k + I - l)Wk{p) 
and, combining with (|49p . we obtain 

{Hx ■ Wk,i{x)) = -{2 -2k- l)——Fu,i. (54) 



dV{p) ' " " ' 'dV{p) 

We therefore obtain the final answer for the free energy: 

which enables one to calculate all the terms J-k,i except the contribution at /c = 1, / = (which is 
the torus approximation in the Hermitian one-matrix model calculated in [3]) and the second-order 
correction in C (the term Ji),2)- The term J^o,2 was calculated in [5], in the next subsection we present 
the answer. 

To calculate all other terms of the free energy expansion we only need to introduce the new 

jc y{s)ds 

vertex q at which we place the (nonlocal) integral term j>^{^) ^ 2y[() — • '^^ obtain Tg^i in the above 
diagrammatic technique we merely remove the root of the tree and replace the first internal vertex 
(which is always white) by this new (two-valent) vertex, as shown in examples below. 

In the Hermitian matrix model case, it was possible to perform the integration for 2y{s)ds 
through the branch point fia in the vicinity of the ath branch point; here it is no more the case and 
we must consider global integrations. Note, however, that we must introduce nonlocal terms only for 
the very last, outermost, integration; all internal integrations can be performed by taking residues at 
branch points and at the zeros of the polynomial M{p). 

For example, the first correction to the free energy is 



-1 • ^0,1 = ® « 



(56) 



By integration by parts (in which we make substitution at q and q, which cancels the factor 2 in the 
denominator), we obtain 

jq y ^ ^ 
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i.e., we come to the semiclassical Dyson term (the entropy). Other lower-order terms of the free-energy 
expansion may also have a physical interpretation, as we demonstrate on the example of the term -Fo,2- 

For the term T\ \ we have 



-^1,1 = 




where we explicitly crossed out the diagram that is forbidden by the selection rules. So, the actual 
answer comprises five admitted diagrams. 

5.1 Calculating J^o,2 

It was shown in |5] that the diagrammatic representation formally applied to just gives zero. 

We can however make a guess for the actual J^o,2- It was proposed in [5] to interpret it as the 
Polyakov's gravitational anomaly term J / R^R, where R is the curvature (of two-dimensional metric) 
and 1/A is the Green's function for the Laplace operator, which in our case is the logarithm of the 
Prime form. The curvature is expressed through the function y as i? ~ y' /u- That is, we have two 
natural candidates for J-0,2: 

To,2 - [ [ dqdpy^ log E{p,q)y^, 

J J y[q) y{p) 



where E is the Prime form, or 

-^0,2 



j j dqdp\ogy{q)B{p,q) logy(p). 



but neither of these expressions is well defined. The first one develops the logarithmic cut oX, p = q and 
cannot be written in a contour-independent way; moreover, both these expressions are divergent when 
integrating along the support. We therefore must find another representation imitating this term. A 
proper choice turns out to be the one in which we integrate by part only once therefore producing the 
expression of the form 

where the second integral is taken exactly along the eigenvalue support T). 

Variation of the logarithmic term in (I57p can be presented already in the form of the contour 
integral. Referring the reader for details to [5], we present here only the result. 

Introducing M^^ = Q(p)|p=^^ to be the first moments of the total potential V and A(//) to be 
the Vandermonde determinant of the branch points /ia, we have 

-^0,2 =-i,^, (A ^^^^^) ' ^^^^ 

and we indeed have quantum correction term (the second term in 



16 



Acknowledgments 



The author thanks Bertrand Eynard for cohaboration, Sara Pasquetti and Marcos Marino for the useful 
discussion, Gabriel Cardoso for the hospitality during my visit to the Program in Matrix Models in 
Instituto Superior Tecnico, Lisbon, Portugal, and Eyjafjallajokull for ensuring the possibility of this 
long discussion. 

The work is supported by the Russian Foundation for Basic Research (Grants 10-02-01315_a, 
10-01-92104-JF_a and 09-01-12150-ofi_m), by the Grant for Supporting Leading Scientific Schools 
NSh- 795. 2008.1, and by the Program Mathematical Methods for Nonlinear Dynamics. 

References 



[1 



[3] 
[4] 
[5 



[7; 



[9 



L. F. Alday, D. Gaiotto, and Y. Tachikawa, Liouville correlation function from four- dimensional 
gauge theories, Lett. Math. Phys. 91 (2010) 167-197; |arXlvr0906.3219 [hep-th]. 

G. Borot, B. Eynard, S. Majumdar, and C. Nadal, Large deviations of the maximal eigenvalue 
of random matrices, arXiv: 1009: 1945 [math-ph]. 

L. Chekhov, Genus one corrections to multi-cut matrix model solutions, Theor. Math. Phys. 
141 (2004) 1640-1653, hep-th/ 0401089| 

L. Chekhov, B. Eynard, Hermitean matrix model free energy: Feynman graph technique for all 
genera, JHEP 0603:014 (2006). 

L. Chekhov and B. Eynard, Matrix eigenvalue model: Feynman graph technique for all genera, 
JHEP 12(2006)026. 

R. Dijkgraaf and C. Vafa, Toda theories, matrix models, topological strings, and N = 2 gauge 
systems, larXiv:0909.2453l [hep-th]. 

T. Eguchi and K. Maruyoshi, Penner-type matrix model and Seiberg-Witten theory, JHEP 
1002(2010)022; larXiv:0911.4797l [hep-th]. 

B. Eynard, Topological expansion for the 1-Hermitian matrix model correlation functions, JHEP 



0411(2004)031, |hep-th/0407261[ 



A. Mironov and A. Morozov, The power of Nekrasov functions, Phys. Lett. B680, 188 (2009) 
[arXiv:0908,2190J [hep-th]]. 

A. Marshakov, A. Mironov and A. Morozov, On non-conformal limit of the AGT relations 
iarXiv:0909.2052 [hep-th]. 

A. Marshakov, A. Mironov and A. Morozov, Zamolodchikov asymptotic formula and instanton 
expansion in N = 2 SUSY Nf = 2Nc QCD, larXiv:0909.3338l [hep-th]. 



17 



